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ABSTRACT 



New and archival interferometric 12 C0(1— >0) datasets from six nearby galaxies are 
combined with H2 2.122 fim and Ha maps to explore in detail the interstellar medium 
in different star-forming galaxies. We investigate the relation between warm (H2 at 
T ~ 2000 K) and cold (CO at T ~ 50 K) molecular gas from 100 pc to 2 kpc scales. On 
these scales, the ratio of warm-to-cold molecular hydrogen correlates with the 
ratio, a ratio that tracks the star formation activity level. This result also holds for the 
global properties of galaxies from a much larger sample drawn from the literature. The 
trend persists for over three orders of magnitude in the mass ratio, regardless of source 
nuclear activity. 

Subject headings: galaxies: ISM — infrared: galaxies — radio lines: galaxies 



1. INTRODUCTION 

Carbon monoxide is an important probe of the interstellar medium. It is used to map large- 
and small-scale kinematic signatures within galaxies, it traces the sites of possible future star 
formation, and it is a favored proxy for estimating the amount of H2. H2 is the predominant 
interstellar molecular species but it is more difficult to detect than CO. Researchers frequently 
use ratios of various radio transitions of CO to gauge the physical conditions of the interstellar 
medium. Such transitions reflect cool gas at temperatures of order 5 — 20 K for disks and 50 — 75 K 
for nuclear regions (e.g. Young & Scoville 1991; Harris et al. 1991; Hafok & Stutzki 2003). The 
mid-infrared spectroscopic capabilities of the Infrared Space Observatory enabled observations of 
warm H2 gas (T ~ 150 — 700 K) via emission that arises from pure rotation of the molecule 
(e.g. Timmermann et al. 1996; Valentijn & van der Werf 1999; Fuente et al. 1999; Rigopoulou 
et al. 2002). Rigopoulou et al. found that the mass fraction of this warm molecular gas, with 
respect to the total molecular g ELS clS traced by 12 CO(1^0), is 1-10% for starbursts and 2-35% for 
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Seyferts (see also Armus et al. 2004; Rodriguez- Fernandez et al. 2004). The fraction of molecular 
material at high temperatures (T ~ 2000 K) indicates the amount of molecular gas either currently 
participating in star formation or experiencing its immediate after effects such as shocks (e.g. 
Scoville et al. 1982; Shull h Beckwith 1982). A useful probe of molecular gas at high excitation is 
the ro-vibrational H 2 2.122 fxva. line. H 2 2.122 /j,m emission derives from collisional excitation by 
shocks or radiative excitation in ultraviolet radiation fields (e.g. Goldader et al. 1997), or possibly 
from reprocessed X-ray photons (Maloney, Hollenbach, & Tielens 1996). In this work we combine 
12 CO(1^0) and H2 2.122 fim data from six nearby galaxies to estimate the fraction of molecular 
material at relatively high excitation. 

The sample is drawn from the large collection of galaxies studied in the ISO Key Project 
on the Interstellar Medium in Normal Galaxies (Helou et al. 1996; Dale et al. 2000). The data 
obtained for the Key Project sample include mid-infrared imaging and spectrophotometry, far- 
infrared spectroscopy, and Ha imaging. As a follow-up to the Key Project we have recently obtained 
5'/7x 10" H2 2.122 pLm. spectral data cubes and arcminutes-scale 12 CO(1^0) datasets for six galaxies. 
The properties of these particular galaxies are fully described in Dale et al. (2004) along with the 
Ha and near-infrared molecular hydrogen data obtained for that work. The targets for this study 
span a range of environments from H II regions to normal galaxy nuclei and perhaps one heavily 
obscured AGN. Our maps of H 2 2.122 fim emission sample sub-kpc to kpc sizescales, so coupling 
those data with high resolution interferometric 12 CO(1^0) maps enables an exploration of hot 
and cold molecular gas on a local scale. These small sizescale data are then compared to the 
disk-averaged properties of galaxies drawn from the literature. 

2. 12 CO(1^0) OBSERVATIONS AND DATA REDUCTION 

The 12 CO(1^0) observations for NGC 1569, NGC 6946, and UGC 2855 are described in 
Taylor et al. (1999), Heifer et al. (2003), and Hiittemeister et al. (1999), respectively. The 
new data obtained for this work are the 12 CO(1^0) maps for NGC 2388, NGC 4418, NGC 7771 
described below. 



2.1. Owens Valley Radio Observatory CO Observations and Data Calibration 

12 CO(1^0) observations were obtained over the period 2003 February through 2003 October. 
The array of six 10.4 m telescopes was used in its "Compact" configuration for observing one 
track each for NGC 2388 and NGC 7771; NGC 4418 was observed in the "Low" and "Equatorial" 
arrangements for one track apiece. The average single-sideband temperature ranged from 304 to 
378 K for the set of four tracks. Spectral resolution was provided by a digital correlator configured 
with four bands, each with 32 x 4 MHz channels, yielding a total velocity coverage 5v ~ 1200 km s _1 
at the redshifted 12 CO(1^0) line. 
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Nearby radio-loud QSOs were used for gain and passband calibration. The data were calibrated 
using the standard Owens Valley array MMA program (Scoville et al. 1993) and mapped using 
the BIMA MIRIAD package (Sault, Teuben, & Wright 1995). The beam sizes for the different 
observations are listed in Table 1 along with the rms noise level in the source-free channels of 
the cleaned maps. The planets Uranus and Neptune were used as flux calibrators; the absolute 
12 CO(1^0) fluxes are uncertain at the 20% level. 

3. RESULTS 
3.1. Morphology 

Contours of the 12 CO(1^0) emission are overlaid on Ha maps in Figures 1-3. Outlines of 
the regions mapped in H2 are included for reference. The spatial profiles of Ha and CO emission 
for NGC 2388 and NGC 4418 are compact and centered on the nuclei. In UGC 2855, NGC 6946, 
and NGC 7771 the 12 CO(1^0) emission generally tracks the Ha distribution whereas NGC 1569 
shows significant spatial differences in the two lines. Additional details for the latter four extended 
sources are described below. 

3.1.1. UGC 2855 

The 12 CO(1^0) emission for UGC 2855 follows the spiral arm and bar patterns seen in Ha. 
The molecular bar is ~8 kpc long. As described in Hiittemeister, Aalto, & Wall (1999), the radio 
continuum and 12 CO(1^0) both peak within an arcsecond of 03 h 48 m 20.8 s +70°07'58". The Ha 
and H2 emission also peak within an arcsecond of this position. The northwestern part of the bar 
is a factor of 1.2-1.5 brighter in 12 CO(1^0) than the southeastern portion (Hiittemeister, Aalto, 
Sz Wall). The apparent linear cutoff along the northern edge is an artifact of the 12 CO(1^0) map 
obtained by Hiittemeister, Aalto, Sz Wall. 

3.1.2. NGC 1569 

Taylor et al. (1999) have previously compared 12 CO(1^0) emission with Ha for NGC 1569. 
We have obtained H2 data for H II region 2 of Waller (1991), which also overlaps with region C of 
Greve et al. (1996) and GMC 3 of Taylor et al. (1999). Neither of the super starclusters for which 
NGC 1569 is famous is covered by the H2 map, though super starcluster A is near the eastern 
edge of the H2 map. The Ha image exhibits fairly striking features due to the winds from the 
starclusters. The H2 emission is strongest in the southwestern quadrant of the region mapped in 
H2 (Dale et al. 2004), a region that also shows significant 12 CO(1^0) emission. 
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3.1.3. NGC 6946 

Sheth et al. (2002) have previously compared the 12 CO(1^0) emission with Ha for NGC 6946. 
The 12 CO (1— >0) emission in the center of NGC 6946 displays a north-south structure and generally 
follows what is seen at Ha, though as pointed out by Sheth and collaborators, the Ha emission 
appears on the leading side of the 12 CO(1^0) emission. 

3.1.4. NGC 7771 

The linear structure to the 12 CO(1^0) emission for NGC 7771 follows the general inclined- 
disk trend evidenced by Ha (and broadband optical and infrared images). As for UGC 2855, the 
molecular disk is asymmetric with the northeastern side of the nucleus 50% stronger in 12 CO(1^0) 
emission than the southwestern. Curiously, the molecular disk emission in the northeast is aligned 
with the Ha emission whereas the two distributions are largely misaligned for the southwestern 
portion of the disk. The Ha morphology essentially follows that at 7 fim (Dale et al. 2000), so the 
northeast-southwest differences cannot be explained by dust extinction. 

3.2. Flux 

The 12 CO(1^0) fluxes integrated over the interferometer maps are listed in Table 2 along with 
fluxes extracted over the field-of-view of the Palomar Integral Field Spectrometer H2 pointings. 1 
In comparison with the global fluxes obtained from single-dish observations that are listed in the 
rightmost column of Table 2, the integrated fluxes from the arrays are systematically lower. Part 
of this discrepancy can be attributed to the relatively small 12 CO(1^0) map sizes (for NGC 6946 
and UGC 2855). Some of the extended flux is also missed by the limited uv coverage inherent to 
interferometric measurements; the minimum dish spacing of the OVRO "compact" configuration, 
for example, resolves out structures larger than ~ 35" (see Heifer et al. 2003). Detailed inves- 
tigative work by Heifer et al. shows that, on 10" scales, interferometers capture 80-90% of the 
total emission. The recovery fraction approaches 100% on these scales if the CO is dominated by 
a central concentration, similar to the situation for all our targets save for NGC 1569. 

3.3. Warm and Cold Molecular Hydrogen Masses 

Rigopoulou et al. (2002) find that the percentage of "warm" molecular gas as traced by mid- 
infrared H2 rotational lines (T ~ 150 K), with respect to the bulk of the molecular gas that is 
"cold" and traced by the millimeter 12 CO(l— >0) transition, is 1-10% by mass for starbursts, and 



iNote that the CO beams are ~30% larger than the PIFS fields-of-view for NGC 2388 and NGC 7771. 
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2-35% for Seyferts. The H2 mass at even higher temperatures (T ~ 2000 K) can be computed from 
the extinction-corrected 2.122 ^urn ro-vibrational H2 near-infrared line flux: 

M(H 2 ) wa rm = 2m p n H2 V H2 = 7 - A — 7- - 5.08 M — i — — 1) 

_A,=i,j= 3 4,=i^0S(i)^ VlO- 16 W m 2 J \MpcJ 

where we assume the Orion vibrational temperature T vi b = 2000 K and thus a population frac- 
tion f v= i J=3 = 0.0122 (Scoville et al. 1982), and a 1— >0 S(l) transition rate of A/=i->oS(i) = 
3.47- 10~ 7 s- 1 (Turner, Kirby-Docken, & Dal garno 1977). The "cold" molecular hydrogen mass can 
be inferred from the integrated 12 CO(1^0) flux according to 

where «g is the Galactic 12 CO(1^0)-to-H 2 conversion factor (Wilson 1995). The warm/cold H2 
mass ratios range between 10~ 7 to 10~ 5 for the small fields mapped in H 2 2.122 /j,m. 



4. DISCUSSION AND SUMMARY 

The luminous infrared galaxy NGC 6240 has the largest known H 2 2.122 fim line luminosity, 
and the largest ratio of H 2 line-to-bolometric luminosity ratio (Lh 2 2.122 ~ 7.5 x 10 7 L , which is 
0.3% of the bolometric luminosity; see Draine k, Woods 1990). The upper end of the range we find 
for our sample's local warm/cold molecular mass ratio (10~ 7 to 10~ 5 ) is consistent with the value 
for NGC 6240 (van der Werf et al. 1993; Young et al. 1995). To expand upon our result and place 
it within a larger context, we have extracted a large pool of H 2 2.122 ^m and 12 CO(1^0) data 
from the literature. A large, uniform database of global 12 CO(1^0) fluxes for nearby galaxies is 
provided by Young et al. (1995). H 2 2.122 fim data for at least 32 of these galaxies can be found 
in the literature, and for those galaxies with multiple measurements we use the fluxes from the 
largest apertures (Hall et al. 1981; Rieke et al. 1985; Heckman et al. 1986; Kawara et al. 1987; 
Fischer et al. 1987; Moorwood & Oliva 1988; Puxley et al. 1988; Lester et al. 1990; Forbes et 
al. 1993; van der Werf et al. 1993; Calzetti 1997; Larkin et al. 1998; Roussel et al. 2003). Since 
the H 2 2.122 ^m maps from the literature typically do not span areas as large as those covered 
by the Young et al. CO maps, we have estimated the portion of the CO flux that derives from 
within the H 2 2.122 apertures. We use an exponential model for the CO emission along with 
the scalelengths and effective CO radii suggested in Young et al. The data are plotted in Figure 4 
as a function of global IRAS far-infrared color. The data are corrected for Galactic extinction and 
a uniform a/ac = 1 is assumed. 

The range in "global" warm /cold H 2 mass ratio drawn from the literature spans a slightly larger 
range, extending down to 10~ 7 - 5 . In addition, the mass ratio tracks the far-infrared color j^^^y- 
Note also that the trend holds regardless of nuclear activity. Correcting for internal extinction would 
likely steepen the slope of the trend, as sources with warmer far-infrared colors typically exhibit 
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larger extinctions (e.g. Hattori et al. 2004). However, compared to the Ay < 3 mag typically 
found for starburst nuclei (Ho, Filippenko, & Sargent 1997; Dale et al. 2004), the reddening vector 
displayed in Figure 4 shows that a large internal extinction is necessary to significantly shift the 
data points. Another caveat is that the CO-to-H2 conversion factor is likely to differ from galaxy to 
galaxy. Some authors suggest the CO-to-H 2 conversion factor clearly depends on metallicity (e.g. 
Wilson 1995; Arimoto, Sofue, Tsujimoto et al. 1996; Israel 1997) whereas others claim negligible 
variation above a metallicity of 12+log[O/H]~8.0 (e.g. Rosolowsky et al. 2003; Bolatto et al. 
2003). Metallicities were found for half of the galaxies in the literature, and using the CO-to-H2 
prescription described in (Wilson 1995) we find only a 10% increase (0.04 dex) in the average 
warm/cold ratio. If internal extinction and the CO-to-H2 conversion factor a were properly gauged 
for each galaxy the data points would differ, but it is unlikely that the main trend would be washed 
out. 

It is worth questioning whether this trend with jjftxfj^ is just a global effect: Are the vari- 
ations from the overall trend more significant on a local scale? Does the observed trend simply 
reflect smoothing on a global scale of more variable local phenomena? The 12 CO(1^0) interfero- 
metric maps allow us to sample the cold molecular hydrogen on a scale more similar to the typical 
extent of the H2 2.122 fim emission. The local far-infrared colors are estimated using ^(15^^ 

from the higher resolution ISO mid-infrared imaging data and the ^^^^ -to- j-^^qq^^ mapping 
for star-forming galaxies (Dale et al. 2001; Dale & Helou 2002). The small data points displayed in 
Figure 4 represent the "local" ratios for the six galaxies in Table 1 — these data stem from just the 
Palomar Integral Field Spectrometer H2 fields-of-view. Except for NGC 4418, which is listed as a 
Seyfert 2 galaxy in NED, the local ratios are indicated by stars to indicate that they derive from 
star-forming regions. The data for these six small-scale (sub-kpc to kpc) regions are consistent 
with the galaxy global data. The relation may break down on the parsec or tens of parsec scale 
(e.g. Scoville et al. (1982), where the mass ratio is 10~ 3 — 10~ 2 for the core of the Orion molecular 
cloud), but higher resolution data would be required to probe to this level in extragalactic regions. 

The gas mass ratio gauges the molecular gas phase of the interstellar medium whereas the far- 
infrared color is a probe of the dust emission. The two measures are coupled, likely due to the effects 
of star formation on the interstellar medium. The warm/cold molecular gas mass ratio is a measure 
of star formation activity via shocks and perhaps UV excitation (Draine & Woods 1990; van der 
Werf et al. 1993), while the far- infrared color is a measure of star formation activity via dust grain 
temperature (Hattori et al. 2004; Xilouris et al. 2004; Dale et al. in preparation). In either case, 
the same stars are ultimately responsible for heating the warm molecular gas and interstellar dust. 
A locally-based coupling could arise from molecule formation on dust or molecular gas excitation 
via the photoelectric effect from grains. On the other hand, a more global description would simply 
link radiation leaking out from a younger stellar population embedded within or near a molecular 
cloud complex and heating interstellar dust distributed on kpc scales. The photoelectric basis for 
local coupling is unlikely since photoelectrons primarily arise from the smallest grains, whereas the 
far- infrared color is dominated by emission from large grains (e.g. Hollenbach & Tielens 1999). 
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Moreover, the far-infrared color is more of a large-scale measure since it reflects the balance of 
emission from diffuse cirrus and H Il/photodissociation regions (e.g. Helou 1986). In other words, 
a global scaling is likely. 

We have obtained near-infrared H2 and interferometric 12 CO(1^0) maps for six nearby galax- 
ies, data that trace different phases of interstellar molecular gas. We find that the mass ratio of 
warm and cold molecular gas is of order 10~ 7 to 10~ 5 , indicative of the amount of molecular ma- 
terial that is experiencing the effects of active star formation. After supplementing these sub-kpc 
to kpc scale data with integrated (global) data from the literature, we find a broad trend between 
the molecular gas mass ratio and the far-infrared color. Both ratios are conceptually linked to star 
formation, though it is unlikely that they are directly physically coupled. It is more likely that the 
link relates to radiation emerging from young stars within or near molecular clouds and heating 
dust grains distributed on kpc scales. 

We are thankful to C. Taylor for the use of his NGC 1569 CO map. C. Kobulnicky, T. Jarrett, 
and E. Murphy provided useful comments and input. The Owens Valley millimeter array is sup- 
ported by NSF grant AST 99-81546. This research has made use of the NASA/IPAC Extragalactic 
Database which is operated by JPL/Caltech, under contract with NASA. 
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o US 

Fig. 1. — Integrated 12 CO(1^0) emission contours overlaid on continuum-subtracted Ha images 
for UGC 2855 and NGC 1569. The contours for UGC 2855 start at 3 Jy beam" 1 km s' 1 and are 
spaced by the same value; they start at 0.2 Jy beam -1 km s _1 and are spaced by the same value for 
NGC 1569. The ellipses represent the beam sizes for the CO observations. Rectangular outlines of 
the regions mapped in H2 2.122 /im are included for reference; grayscale images of the H2 maps are 
inset in the lower righthand corner of each image. The 12 CO(1^0) and Ha maps for NGC 1569 
are from Taylor et al. (1999) and Martin, Kobulnicky, & Heckman et al. (2002), respectively. The 
12 CO (l->0) map for UGC 2855 is from Hiittemeister, Aalto, & Wall (1999). 
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Fig. 3.— Similar to Figure 1 for NGC 6946 and NGC 7771. The contours start at 
15 Jy beam -1 km s _1 for NGC 6946 and at 10 Jy beam™ 1 km s™ 1 for NGC 7771, respectively, and 
are spaced by the same values. The 12 CO(1^0) and Ha maps for NGC 6946 are from Heifer et al. 
(2003) and Kennicutt et al. (2003), respectively. 
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a OVRO-based 12 CO (1 -*•()) fluxes are uncertain at the 20% level. 
b FCRAO data (Young et al. 1995). 

c Mass ratios from the PIFS regions are uncertain at the 50% level, driven by the large absolute uncertainty in the 
H 2 2.122 fluxes (see Dale et al. 2004). 

^Infrared colors from the PIFS regions are uncertain at the 30% level, based on a combination of uncertainties in the 
infrared fluxes and models (see § 4 and Dale et al. 2000). 
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Fig. 4. — The ratio of warm-to-cold molecular gas mass as a function of far-infrared color. 



